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We determine the spreading resistance for a sheet of homogeneous ma- 
terial of uniform thickness with a disk contact source on one side and with 
current collected (i) over the entire back plane, and (ii) at a corresponding 
disk on the back plane. A constant driving potential is assumed over the 
source resulting in mixed boundary conditions in the plane of the source. 
For each case, a closed form integral solution is derived and then numeri- 
cally integrated for a range of geometric ratios resulting in a universal 
spreading resistance curve. The results are used to evaluate the spreading 
resistances encountered in a typical (1 Q cm) semiconductor material. 

I. INTRODUCTION 

The resistance associated with the nonparallel current flow between 
a spacially separated source and sink is referred to as spreading 
resistance. Calculation of spreading resistance is often required in the 
analytical treatment of semiconductor devices. In particular, elec- 
trical current flow in a slice of silicon between a surface contact and 
a back-plane contact involves the calculation of the ohmic spreading 
resistance. The heat flow between an active transistor or integrated 
circuit and an external heat sink involves a calculation of the thermal 
spreading resistance in the device carrier. Also, for a given structure, 
the capacitance including fringing is directly related to the conduct- 
ance including spreading. 

Two cases are considered in this paper. The spreading resistance is 
determined for an infinite sheet of homogeneous material of uniform 
thickness with a disk contact source on one side and with current 
collected (i) at a completely metallized back plane, and (ii) at a 
corresponding disk on the back plane. D. P. Kennedy analyzed 
these cases for finite cylindrical volumes but with nonmixed boundary 
conditions. 1 He assumed a constant flux over the source region and 
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defined spreading resistance in terms of the maximum temperature (or 
potential) on the disk. In this paper, we assume a constant driving 
potential over the source. This assumption results in a mixed boundary- 
condition in the source plane because the potential is specified over 
part of the surface, and the normal derivative is specified over 
the rest of the surface. In the limiting case of a cylinder of infinite 
radius, there appears to be some difference between Kennedy's results 
and those presented herein. This appears to be associated with Ken- 
nedy's definition of spreading resistance in terms of the maximum 
disk potential. However, the comparison is based on the extrapolation 
of curves given by Kennedy. 

A. Gray, et al., considered certain special cases of the infinite region 
problem with mixed boundary conditions in the source plane. 2 How- 
ever, they imposed one of two constraints. Either (i) there was 
sufficient separation between source and sink so the flux distribution 
at the disk could be taken equal to the limiting half-space case or 
(u) the disk was small enough to be considered a point source. 
In essence, these constraints again reduce the problem to one with a 
nonmixed boundary condition in the source plane. Neither constraint 
is imposed in this paper. 

In this paper, for each case, a closed form integral solution is 
derived which is numerically integrated for a range of geometric 
ratios resulting in a universal spreading resistance curve. For both 
cases, the curves approach the half-space limit for thick sheets, and 
the curves asymptotically approach the nonfringing limit for very 
thin sheets. ' Finally, the results are used to evaluate the spreading 
resistances typical of those encountered in semiconductor technology. 

II. ALGEBRAIC SOLUTION 

The geometry and system of coordinates are given in Fig. 1. Solving 
La Place's equation in cylindrical coordinates when the potential <f> is 
independent of 4 gives 3 

<f>(p, z) = I i(k)(Ci cosh kz + C 2 sinh kz)J (kp) dk, (1) 

•'0 

where f(k) must be determined by the boundary conditions at z = 0. 
For nonmixed boundary conditions at 2 = 0, f(k) can be found by 
inverting the Hankel transform. 1 In the present case, however, the 
mixed boundary condition must be imposed at the 2 = surface. Thus, 
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d0 
dz 



V = T Kk)dJo(kp) dk, fg p £ a; 

= I /(fc)fc(Ci sinh kz + C 2 cosh kz)J (kp) dk , 

«-0 Jo 

This results in a dual set of integral equations. 

V = T j(k)C,J {kp) dk, ^ p £ a; 

•'0 

= [ j(k)C 2 J (kp) dk, p^a. 
Jo 



(2) 



= 0. 



(3) 



J. D. Jackson observes that the solution of this set of equations is 3 

W -S§*fe (4) 

Thus 

^(p ( z) = / — 7 — (cosh kz + Trsinh kzJJ (kp) dk. (5) 

2.1 Case ./ — Back Plane Grounded 

In the first case considered, the second electrode is a completely 
metallized or grounded back plane. The boundary conditions are 
given in Fig. 2. This situation is representative of heat flow from 
an integrated circuit through an insulated header to a can acting as 



DISK AT CONSTANT 
POTENTIAL, V 




Fig. 1 — Physical geometry and coordinate system for calculating spreading 
resistance. 
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--HALF-SPACE LIMIT 
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Fig. 2 — Spreading resistance i?i between disk and grounded back plane. 



a thermal radiator. 



At z = w, 4>{p, z) L_ I(1 = 0, 



and thus 



/ ' - '" '"* J (kp)[ cosh kw + ^rsinh kw ) dk = 0. 



Since this must hold for all w, 



C 

jf sinh kw = — cosh kw. 



Thus, 



<f}(p, z) = - — / —7— JJkp)(cosh kz — coth kw sinh fa?) dk. 
t J a ka 



Now, 



./ = aE = — o-V0 at z = 0, 
1=1 I Jpdp dd, 

Jo Jo 



(6) 

(7) 

(8) 

(9) 

(10) 
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and 



V4> 



dz 



= -— [ «m ka,J (k P ) coth kwdk. (11) 

Combining the last throe equations and carrying out the indicated 
integrations, yields 



I = IV aa 



f 

•'o 



co th kw sin kaJ^ka) 
k 



dk. 



Thus 



and finally, 



1 l i 



* coth fcig sin kaJi(k a) 
k 



dk 



(12) 



(13) 



#,<ra = 



4 f s -^ ./,(.,) coth (=)* 

«/ 1) .( \ (I / 



(14) 



2.2 Case £ — Z^'s/i 1 on Back Plane Grounded 

The case where the sink consists of a disk of radius a coaxial with 
the source can be derived by image theory from Case 1. A drawing 
of this configuration is shown in Fig. 3. The values of spreading resist- 
ance in Case 2 can he derived from Case 1 by setting 



7?,(7fl = 



Rooa 



(15) 



and 



2a 
w 



The resulting spreading resistance equation \> 

1 



R.,aa = 



■jsp '.«-*©* 



(16) 



(17) 



III. NUMERICAL EVALUATION 



3.i Numerical Integration — Universal Curves 

The integrals in equations (14) and (17) can be evaluated numer- 
ically for various a/ic ratios if the upper limit is finite and the 
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Fig. 3 — Spreading resistance R a between disk and grounded disk on back plane. 



integrands behave well at the lower limit. Investigation of the be- 
havior of the integrands over the entire range of integration, indicated 
the range could be split into three ranges: to 0.1, 0.1 to 300, and 
300 to oo. For the 300 to oo range, and a/w > 1, the integrals for both 
cases reduce to 



I 



sin x 



300 # 



J x (x) dx = 1 3 . 



But 



or' 



^oo • i*30Q 

T I S1TiX T t \ A I SmX T 1 \ J 

h = I «/iW dx — I J\\x) dx 

JQ X J § X 

/»300 • 

h = 1 - / S -^^ Jx{x) dxml-It. 

Jq x 



(18) 



(19) 



(20) 



7 4 was evaluated numerically on a digital computer and found to be 
0.96439. Thus Z 3 = 0.03561 and is independent of a/w for a/w ^ 1. 

For the to 0.1 range, the integrands were replaced by their small 
argument approximations. For x <C 1, 



SPREADING RESISTANCE 781 



J.W-jg-f. (2D 



sin x 
x 



1, (22) 



eofcHSUJL, ( 2 3) 

a wx 



and, 



., tux 2a , . x 

cothr > (24) 

2a iwx 

It follows that the small argument approximations for the integrands 
are a/2iu, and a/w for Cases 1 and 2, respectively, and the values of 
the integrals over the range to 0.1 are a/20w and a/10w f respectively. 
This is a good approximation because the integrands are very flat 
functions of x near x = 0. 
Thus for numerical purposes, 

R^ = -r =~. " 7-x =i (25) 



fe + r s ^^ )cotii ( 2 ?) <fa+o - 03561 ] 



and 



R * aa = "f 7^~. 7 — \ =1 ' ( 26 ) 

2 Li^ + /,, s f '•« coth (1) * + oo356i J 

These results hold for a/w i^ 1. For a/w = 0, the half-space limiting 
case for finite a and infinite w can be used to find values for Ri and 
R 2 . The resultant values are I and £ for cases 1 and 2, respectively. 
This result follows directly from Jackson's work. 3 

Using the half-space limits and carrying out the numerical integra- 
tions for various values of a/w results in the universal spreading 
resistance curves given in Figs. 2 and 3. 

3.2 Asymptotic Limits for Cases 1 and 2 

As the ratio a/w becomes larger, the components of R i and R 2 due to 
fringing become progressively smaller. Neglecting fringing, the resistance 
for either case is 

R - -£? (27) 
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or 



Roa = — • 

TTd 



(28) 



The results given in this section can be presented more clearly in 
terms of conductances because the fringing and nonfringing components 
are in parallel and RT 1 = G is a linear function of a/w. Let the con- 
ductances corresponding to Ri and R 2 be designated Gi and G 2 , re- 
spectively. For a given a and w, there should be more fringing in Case 1 
than in Case 2, and no fringing for R. Thus for any a/w, 



G t > G 2 > G. 



(29) 



Also Gi and G 2 should asymptotically approach Gasa function of a/w. 
These observations are supported by the curves of G x , G 2 , and G 
versus a/w given in Fig. 4. The differences of fringing components for 
Cases 1 and 2 are given in Fig. 5. Gi , G 2 and G are within one percent 
for a/w ^ 10. 

3.3 Typical Example 

The results given above have been used to evaluate spreading resis- 
tances typical of those encountered in semiconductor technology. A 
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Fig. 4 — Spreading conductances with and without fringing between disk and 
grounded back plane and between disk and grounded disk on back plane. 
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Fig. 5 — Spreading conductance fringing components. 
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Fig. 6 — Spreading resistances .fti and R a for 1 ft cm material. 
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1 fi cm slice of silicon with circular contacts was considered. Figure 6 
gives the calculated #! and R 2 for such a configuration. 

IV. SUMMARY AND CONCLUSIONS 

An investigation was made to determine the spreading resistance in 
a sheet of homogeneous material of uniform thickness with a disk 
contact source on one side and with current collected (i) at the entire 
back plane, and (u) at a corresponding disk on the back plane. These 
cases were analyzed and evaluated exactly by solution of a dual set 
of integral equations. The method represented the boundary conditions 
as they physically exist. In contrast to previous work, a single, uni- 
versal, spreading resistance curve was presented for each case. These 
curves should be useful in designing devices and analyzing materials. 
In particular, the curves could be used to determine conductivity of 
a sheet of semiconductor material if its thickness is known. Also, the 
curves could be used directly in the calculation of total capacitance 
and fringing capacitance by relabeling the ordinates with ea/C instead 
of Raa and C/ta instead of G/ua. 
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